Probing the symmetry energy from the nuclear isoscahng 
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Using different parameterizations of tlie nuclear mass formula, we study the sensitivity of the 
isoscaling parameters to the mass formula employed in grand-canonical calculations. Previous works 
on isoscaling have suggested that the symmetry energy implied in such calculations is anomalously 
smaller than that suggested by fits to nuclear masses. We show that surface corrections to the sym- 
metry energy naturally broadens the isotopic distribution thus allowing for values of the symmetry 
energy which more closely match those obtained from nuclear masses. 

PACS numbers: 25.70.Pq, 24.60.k 



I. INTRODUCTION 

Understanding multi-fragmentation, the breakup of 
nuclei into many fragments, has challenged theorists and 
experimentalists for almost three decades [l|, y, 0, 0, Q ■ 
This phenomenon can be observed in both relativistic 
light-ion collisions and in central nucleus- nucleus colli- 
sions, even though the role of reaction dynamics in these 
tvifo domains can be quite different. The consequential 
variation in the mixture of statistical and time-dependent 
phenomena in multi-fragmentation complicates the ap- 
plication of the traditional theoretical tools. Nonethe- 
less, statistical techniques have been successfully applied 
to calculate the relative production of various fragments 
0j 0, S H 01' G^S'n. though the fundamental statistical 
parameters such as temperature or density are not ex- 
ternal constraints, but must be inferred from data (see 
0, HSQ and references therein). 

It is therefore natural to apply statistical methods to 
fragment is otop e distributions and to their isoscaling be- 
havior [l3, [ll|. The latter term refers to the scaling 
properties observed in Refs. [13, [ill, that govern the ex- 
perimentally observed ratios between the yields Yi {N, Z) 
associated with an emitted fragment of atomic and neu- 
tron numbers Z and iV, respectively, produced in two 
similar reactions, labeled i = 1,2: 



R2i{N,Z) = Y2{N,Z)/Yi{N,Z)o^G^X^{aN + l3Z) , (1) 

where a, and /3 are the isoscaling parameters. The 
isoscaling of fragment isotopic yields has also been 
observed in deep inelastic collisions, evaporation pro- 
cesses, as well as in nuclear multifragmentation reactions 
[111 . \VA [l3j . It has been exploited by many authors to 
obtain information on observables related to the multi- 
fragment emission [ll|, [iJ, [la, [la, [l3l- In particular, a 
link between the symmetry energy coefficient Cgym and 
the isoscaling parameter a has been suggested in Refs. 
0,[ll, 
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If a reaction allows straight-forward application of statis- 
tical approaches, it may be possible to obtain empirical 
constraints on Cgym, at different values of the breakup 
density and temperature T, from the knowledge of the 
decaying sources {Ai,Zi) and {A2, Z2) or, the proton 
fractions, Zi/Ai and Z2/A2, together with a and T. To 
explore this possibility, many studies al ong this line have 
been carried out, which rely directly [l3| or indirectly 
[iJ, [la, [H, [20, ^M on the close connection between the 
symmetry energy and the isoscaling parameter a. Such 
studies require explicit assumptions about the form for 
the symmetry free energy and about how the symmetry 
free energy depends on temperature. 

In this work, we investigate these assumptions. We ex- 
plore the sensitivity of a to the mass formula employed 
in statistical calculations. In addition to the standard 
parametrization of the binding energy [l4| , we also con- 
sider mass formulae in which the surface effects in the 
symmetry energy are included |22l . [23| . We examine the 
rigor of present extractions of the temperature depen- 
dence of the symmetry energy from isoscaling analyses 
and whether such analyses presently provide unambigu- 
ous information on this temperature dependence. In sect. 
im we present the statistical treatment employed in this 
work, as well as the mass formulae we use. The sensitiv- 
ity of the isoscaling parameters to the binding energy is 
then investigated in sect. IIIII Our main conclusions are 
delineated in sect. IIVI 



II. THE MODEL 

The Statistical Multifragmentation Model (SMM) ^, 
[25l . [2q is based on a scenario in which a hot source with 
mass and atomic numbers Aq and Zq, respectively, at 



temperature T, expands to a density p ~ po/(l + x) 
and undergoes a prompt breakup. The normal nuclear 
density is denoted by po and 1/x corresponds to the ratio 
between the breakup volume and that available to the 
fragments. The properties of the system are calculated at 
the freeze-out stage, where it is explicitly assumed that, 
except for the Coulomb repulsion, there is no interaction 
among the fragments. Some interactions of the fragments 
with the surrounding light particle gas are incorporated 
via the temperature dependence of the nuclear surface 
energy. Within the SMM approach, the total Hclmholtz 
free energy of the system is approximated by 
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where rii denotes the multiplicity of species i, fi and 
jtrans j-espectivcly stand for the internal and translational 
contributions of a fragment of species i to the free en- 
ergy. The last term in the above expression corresponds 
to the homogeneously charged sphere which, with the 
other terms contained in fi, give the Wigner-Seitz ap- 
proximation to the Coulomb interaction among the frag- 
ments [2J, [231 . The coefficient Cc is taken from a mass 
formula (see below). The species multiplicities are con- 
strained by the properties of the decaying source: 



The average fragment distribution may be obtained 
by minimizing the total free energy of the system with 
respect to {rifc}, subject to the constraints dH)-©: 
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where 



T = F + flB 
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and p,Q and p,B are the Lagrange multipliers, which turn 
out to be the charge and baryon chemical potentials, re- 
spectively. Thus, the above relations lead to 
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Upon using Eqs. (|4])-([5|), one finally arrives at 



p-sAk - ^.qZk ■ (9) 
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A. The grand canonical approach 

For clarity, we use the grand-canonical version of SMM, 
from which analytical expressions for the average frag- 
ment multiplicities can be easily obtained. 

By construction, the terms involving fi (T, V) have no 
explicit or implicit dependence on rii, since it has con- 
tributions only from the binding energy Bi, the inter- 
nal excitation energy of the fragment and the remaining 
Wigner-Seitz Coulomb terms of Eq. ^. Thus, at this 
point, one only needs the explicit expression for f^'^^'^^ 
which; through the Stirling's formula, log 72! « nlogn— n, 
can be written as 
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where A^ = y27r^?/Wr, m is the nuclcon mass, gi de- 
notes the spin degeneracy factor, and Vf stands for the 
free volume. 
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It should be noticed that the last two terms in the 
above expression arise from the homogeneous term of the 
Wigner-Seitz approximation. Since the chemical poten- 
tials can be redefined to absorb these terms, they do not 
play a role in the determination of the particle multiplic- 
ity. However, they modify the relation between a and 
the mass formula, inasmuch as the difference between 
the chemical potentials in reactions 1 and 2 is the rele- 
vant quantity to the calculation [hI, [ij]- Nevertheless, 
since the considered sources are usually similar, the sym- 
metry energy still gives the main contribution to a and, 
therefore, this correction is neglected. 

The internal free energy of the species i is given by 
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where Bi denotes the binding energy of the fragment, 
eo = 16 MeV, Po = 18 MeV and T^ = 18 MeV. The 



last term in the above relation is associated with the 
Wigner-Seitz correction, whereas the Coulomb self en- 
ergy of the fragment is already included in Bi . The form 
of the suppression of the surface energy as T — > Tc was 
determined in Ref. |24i, b y matching the microscopic re- 
sults obtained in Ref. [28[. Very light fragments, A < 5, 
are considered as point particles and therefore contribute 
to fi only through their Coulomb and binding energies, 
where the latter are taken from empirical data, as well 
as their spin degeneracy factors gi. The exception is the 
alpha particle, in which case, one suppresses the surface 
contribution to the excited states. For all heavier nuclei, 
A > 5, Bi is computed through a given mass formula 
and it is assumed that gi = 1 as the contribution to the 
statistical weight is, to some extent, taken into account 
by the excitation part of the free energy. 

The chemical potentials /ig and /i^ are obtained by 
imposing the constraints given by Eqs. ^ and ^. In 
the numerical calculations, the relative error on the mass 
or proton numbers is required to be better than 10^^'' in 
each case. 



B. The binding energy 

The binding energy of nuclei plays an important role 
in the determination of the relative yields produced in 
the breakup of a nuclear system ^^ and therefore it af- 
fects the isoscaling parameters. To investigate this point, 
we consider three different mass formulae based on the 
Liquid Drop Model (LDM) [1^: 



Ba.z = C,A - CsA^I^ - Cc 



where 
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where N = A — Z denotes the neutron number. 

The symmetry energy has contributions from the sur- 
face and bulk terms, which can be easily identified by 
grouping the corresponding factors from the formula: 
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where 



Cs- 



.41/3 



(17) 



The simplest form of this mass formula, in which fc^ = 0, 
Cd = 0, has been used in many statistical calculations 
0,13 ^^'^ ^ill b*^ labeled henceforth as LDMl, whereas 
the complete formula, used in Refs. [23,1301, is denoted 
as LDM2. Since the contribution associated with Cd is 
small, the main difference between the two formulae is in 
the inclusion of surface effects in the symmetry energy in 
LDM2. 

The third mass formula used in this work, LDM3, was 
introduced in Ref. [24I, which also considers surface ef- 
fects in the symmetry energy, but has a different func- 
tional form: 



Ba.z - aA-a,A2/3 
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This form reflects the well-known correlation between the 
bulk symmetry energy and the dependence of the surface 
energy on the isospin asymmetry [3l|, [33, [3J] . Note that 
a* and /3* do not correspond to the isoscaling parameters. 
We chose these symbols so as to keep the notation as 
close as possible to that used in the original work. The 
Coulomb parameter A, in this case, reads 
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and i = v^s corresponds to volume and surface terms, 
respectively. Notice that Cv and Cs might include con- 
tributions associated with the symmetry energy, as we 
discuss below. The following two terms correspond to 
the ordinary Coulomb reduction of the binding energy 
and its correction due to the surface diffuseness, respec- 
tively. In the first two models we consider below, A = 0, 
so the coefficient Cc = ac is constant. This is not the case 
for the last mass formula discussed below. The pairing 
coefficient has the usual meaning: 
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Since this binding energy formula is more complex than 
that given by Eq. (fT2|) . the symmetry energy cannot be 
cast in the same form as Eq. ([T5|) . For large A and small 
asymmetry, it can be approximated by Eq. P^ . where 
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which is identical to the result found above. 

As pointed out in Ref. [24I, the inclusion of the nor- 
mal surface and surface symmetry energies in the mass 
formula should be expected as it simply reflects the de- 
pendence of the nuclear energy functional or equation of 
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FIG. 1: (Color online) Difference between the calculated bind- 
ing energies and the empirical values of Ref. [3J| as a function 
of the mass number for the three mass formulae used in this 
work. 



state on the nuelcar density. The magnitudes or func- 
tional forms of the surface symmetry energy, however, 
depend on the density dependence of the symmetry en- 
ergy, which is poorly constrained. Investigations of the 
temperature dependence of the symmetry energy implic- 
itly assume the temperature dependence to be a result 
of the underl ying density dependence of the symmetry 
energy [ll|, [l^, [l^l ■ To model a temperature dependence 
without considering the surface symmetry energies of nu- 
clei is intrinsically inconsistent. 

The parameters of the three mass formulae have been 
fitted to the empirical binding energies given in Ref. [3J| , 
for A > A. In fig. [1] we show the difference AB between 
the calculated Ba^z and the empirical values. In all the 
cases, one observes a fairly good overall agreement with 
the data, except in the regions where shell effects are 
relevant and cannot be taken into account by these simple 
formulae. It is important to emphasize that AB denotes 
the total difference and that it is not divided by the mass 
number, as this quantity is usually displayed. Larger 
discrepancies are observed for light nuclei, A < 25, with 
the LDMl than for LDM2 or LDM3, because surface 
effects in the symmetry energy, which are more important 
in the light fragments, is disregarded in the LDMl mass 
formula. The other mass formulae, LDM2 and LDM3, 
are both more accurate than LDMl, and LDM2 is slightly 
more accurate than LDM3. 

The best fit parameters arc listed in table HI Although 
paring is switched off in the calculations presented in the 
next section, we include it here for completeness. One 
may notice that the common parameters in LDM2 and 
LDM3 are very similar. The best fit value of a* is very 
close to ttyky, which explains the resemblance between 
the predictions of the two mass formulae. Deviations 
should be expected for light masses, where the factor 



proportional to 1/A^^^ in Eqs. ([17]) and ([20]) becomes 
non-neghgible, because a*^//3* « 44.3 MeV differs ap- 
preciably from Usks « 33.7 MeV. 



III. RESULTS AND DISCUSSION 

The grand canonical model presented above is now ap- 
plied to the description of the breakup of the ^^^Sn and 
the ^^^Sn nuclei at p/po = 1/3. In Fig. [51 we show the 
isoscaling parameter a, obtained using the different mass 
formulae, as a function of the temperature. At T=6 MeV, 
up to 30% reduction of the isoscaling parameter a can be 
obtained by including the surface effects in the symmetry 
energy (LDM2/LDM3) (see Fig.[l]), while preserving the 
good agreement with the empirical binding energies. The 
final values for a are comparable to the data of ref. [ij] , 
which suggest that one may be able to reproduce data by 
adopting a more accurate mass formula without grossly 
changing the symmetry energy values obtained from the 
fitting of the nuclear binding energies. Even so, the re- 
sults reveal that distinct paramcterizations of the surface 
symmetry energy energy do lead to appreciably different 
values of a, in spite of the fact that the mass formulae 
are, to a large extent, equally accurate. So the param- 
cterizations adopted here may still not be the optimal 
ones. 
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FIG. 2: (Color online) Isoscaling parameter a obtained using 
different binding energies in the grand canonical calculations 
as a function of the temperature. 

This differences in the predictions for a for the dif- 
ferent mass models may be understood from the rela- 
tionship between a and the synmictry energy in Eq. ([2|). 
as discussed in Refs. [llliJJj]. At low temperatures, the 
chemical potentials are strongly influenced by the heavy 
fragments, which dominate the mass distribution. Since, 
in this mass region, the coefficient Cgym in LDM2 and 
in LDM3 is nearly the same, as just discussed, Eq. ^ 
predicts that the isoscaling parameter a obtained with 



TABLE I: Parameters of the liquid drop mass formulae used in this work. All the values are given in MeV. 



Label 
LDMl 
LDM2 
LDM3 



15.2692 
15.6658 
15.5586 



as 
16.038 
18.9952 
18.2369 



fflc 

0.68698 
0.72053 
0.70521 



0.0 

1.74859 

na 



Op 
11.277 
10.857 
11.012 



dyKy 

22.3918 

27.7976 

na 



QjsKs 

0.0 

33.7053 

na 



a 

na 

na 

27.1508 



(3* 

na 

na 

16.6412 



these two mass formulae should be very similar. The op- 
posite holds in the case of the LDMl. One should note 
that Eq. ^ has been derived at very low temperatures 
[14| . As the temperature increases, more light fragments 
with A < 50 are produced and the changes in the chemi- 
cal potentials in the -^^^Sn and ^-^"^Sn systems are appre- 
ciably affected by them. Therefore, the binding energy 
differences observed in Fig.[Tl in this lighter mass region, 
become very important. This explains the increase of 
the discrepancy between the values of a obtained using 
LDM2 and LDM3 at high temperatures. Eventually, our 
three curves will converge again at much higher tempera- 
tures (r >> 10 MeV), where only very light particles are 
produced. This latter convergence would occur because 
we use empirical binding energies in this mass region in 
the statistical calculations in all the cases. 

The difference between the values of a obtained us- 
ing LDMl and LDM2/LDM3 is important for interpret- 
ing experimental data because in recent studies the sym- 
metry energy coefficient has been adjusted in order to 
reproduce the experimental values of a [I^l- Our re- 
sults indicate that a large temperature dependence to 
the symmetry energy may not be necessary if a suitable 
parametrization of the binding energies are incorporated 
in the models. Nevertheless, we continue below to discuss 
another issues that may be important to consider if it still 
appears interesting to include an additional temperature 
dependence to the symmetry energy. 

We now examine the consequences of changing Cgym 
in statistical calculations, as has been made in dif- 
ferent studies that focused on the isoscaling analysis 
[13, [H, [H, [11 and on the isotopic distribution [3i,[33. 
The top panel of Fig. [3] displays the parameter a ob- 
tained using LDMl for different values of Csym = ciyky, 
for T — 6 MeV. The full line represents the results ob- 
tained in the approach of [ij, [H, ]1B, ^M ^'^d [3y, l37j . 
case (a), where all the parameters of the mass formula 
are kept fixed, except for fcg. The values for a decrease 
as the Csym decreases, as noted by Refs. [l7l [isl. [l9l. [35| . 
The left side of the lower two panels in Fig. 3 display the 
difference AS between the calculated Ba.z and the em- 
pirical values. Reducing the symmetry energy coefficient 



to Cs- 



10 MeV, similar to the values deduced in Ref. 



[iTJ l , increases the binding of the most neutron-rich frag- 
ments of Mass number 100 by up to 100 MeV without 
changing the binding of the symmetric fragments. If the 
authors of Refs. [13, [H, \lB, ^M ^'^^ correct, why does 
it occur? Without changing the nuclear density, such a 
large effect can only be due to changes in the interaction 
component of the symmetry energy; recent calculations 
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FIG. 3: (Color online) Upper panel: isoscaling parameter a 
obtained using LDMl as a function of Csym. Lower panels: 
AB for different values of Csym- For details, see the text. 



of the temperature dependence for interaction compo- 
nent of the symmetry energy predicts a a reduction in 
Csym of less than 5 MeV [11, [la, [i3l ■ 

Indeed, the authors of [ll,[ia,[H, explain their reduc- 



tion in Csym to be the result of a reduction in nuclear 
density. This plausible explanation for a reduced value 
for Csym does not accurately re pre sent the calculations 
that were performed by Refs. [ll,[i^,[33|, however. While 
it is true that the overall density of the multi- fragmenting 
systems studied by Refs. [l^ \^ [S^ is reduced, accord- 
ing to the descriptions given in these studies and in the 
references they cite, the increased volume is mainly due 
to the separation between the fragments. The fragments 
in their SMM calculations have the bulk, surface and 
Coulomb energies of nuclei at normal density, and the 
excluded volume corrections in their calculations assume 
all of the fragments to be at normal density. 

We note that it is possible to reduce the overbinding for 
the heaviest neutron-rich isotopes by refitting the binding 
energy expression used in the free energy calculation of 
hot nuclei to the measured masses. We have done this 
and show the results as case (b) in the right panels of Fig. 
3. The dashed line in the upper panel of this figure shows 
the dependence of a on Csym to be rather similar to that 
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studies [361 . ISTJ l in order to obtain broader isotopic distri- 
butions. The inclusion of surface symmetry energy terms 
also has a similar effect in this case. We show, in Fig. 21 
the isotopic distribution for some light nuclei, calculated 
for the breakup of the ^^^Sn at T = 6 MeV, using the 
three different mass formulae discussed in this work, with 
the best fit parameters listed in table HI These results 
clearly show that the A dependence of Cgym in LDM2 
and LDM3 leads to wider distributions as Csym becomes 
smaller for light nuclei, due to the influence of the sym- 
metry energy. 



IV. CONCLUDING REMARKS 



FIG. 4: (Color online) Isotopic distribution of fragments pro- 
duced in the decay of the ^^^Sn nucleus at T = 6 MeV. 



for case (a). After refitting, the neutron-rich nuclei with 
A < 30 are overbound by values ranging up to about 
2 MeV/nucleon, indicating a considerable enhancement 
in the attractive potential energy contribution for such 
neutron-rich nuclei. 

The results for both case (a) and case (b) clearly show 
that a is strongly sensitive to Csym, as reported pre- 



viously [17[. For Csym > 22 MeV, there is almost no 



difference whether one changes the other parameters of 
the mass formula, since the symmetry energy becomes 
so important that its contribution to the determination 
of a overwhelms the other terms. On the other hand, 
for smaller Csym values, the other terms of the mass for- 
mula, neglected in the derivation of Eq. ^, become non- 
negligible and therefore the two procedures lead to dis- 
tinct results. 

For Csym below the best fit value, values for AB of 
the order of 2 MeV/nucleon are obtained for light nu- 
clei ^ < 20 at Csym = 10 MeV. We note that such light 
nuclei are the dominant species at high temperatures in 
any SMM calculation, whether for case (a) , (b) or for the 
calculations with LDMI, LDM2 or LDM3 shown earlier. 
Thus, all attempts within the SMM approach to explain 
the values of a reported by [l3, US, U^ [33 involve modi- 
fying the free energies of the light fragments with A < 30. 
The most important theoretically justified modification 
to these light fragment free energies is to include surface 
symmetry energy terms, to reflect the strong influence 
of the nuclear surface on the symmetry energies of light 
nuclei. The option of reducing Csym to low values like 
Csym = 10 MeV, while leading to similar reductions in 
a, does not seem to have a clear theoretical justiflcation. 
More theoretical work is needed to explain why it is rea- 
sonable to expect such a large reduction in Csym and 
to justify the particular form of the reduction chosen in 

Refs. [iiEEISl. 

We note that Csvm has also been reduced in different 



We have analyzed the role played by the mass formula 
used in statistical calculations to determine the isoscal- 
ing parameter a, as well as in the isotopic distribution 
of fragments produced in the breakup of a thermally ex- 
cited nucleus. Our results show that the reduction of the 
symmetry energy values [13, [iM ilB, [Sa, ^3, [231 1 which 
best fit the empirical masses, may not be necessary to 
reproduce the experimental a and the isotopic distribu- 
tions. Rather, the observed experimental reduction of 
the isoscaling parameter [13, [l3, [131 niight be reproduced 
if the surface corrections to the symmetry energy are in- 
cluded in the mass parametrization used in the statistical 
multifragmentation models. Since non-negligible devia- 
tions from the optimal value of Csym lead to large discrep- 
ancies between the theoretical and the empirical binding 
energies and free energies used in the model calculations, 
such changes should not be taken lightly especially if sim- 
ilar effects can be reproduced by the inclusion of surface 
effects in the mass formula. Certainly, some additional 
theoretical work is needed to justify the forms of the mod- 
ifications to the symmetry energies that have been chosen 
in Refs. [l3, [la, [3 • Furthermore, as pointed out in Ref. 
[4ll |. improved treatments to the decay of the primary 
fragments than the usual Weisskopf approach [42, [43| 
also contribute to isotopic distributions with larger tails. 
Finally, we believe that our conclusions will not be af- 
fected by more realistic treatments which may introduce 
a temperature dependence to the symmetry energy since 
the corresponding effects, at a fixed density, are expected 
to be small |3§, SS, EOl • 
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